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I. INTRODUCTION 



In a recent paper ([l])we have studied the behaviour of the constants G, c and A within 
the framework of a flat FRW cosmological model where we consider the effects of a c-variable 
into the curvature tensor. The study was restricted to the perfect fluid case. In such paper 
we arrived to the conclusion that we have only two physical solutions and furthermore both 
solutions follow a power-law time dependence of the physical quantities on the cosmological 
time while in the previous works, where we have not taken into account the effects of 
a c-variable into the curvature tensor, we were capable of obtaining more solutions. We 
would like to emphasize that in such study we obtained as integration condition that both 
"constants", in spite of considering that they vary, verify the relationship G/c 2 = const, for 
all equation of state. Nevertheless we have not been able to find any physical restriction in 
the parameters to determine if the "constants" are growing or decreasing functions on time 
t, since both cases are permissible. We would like to point out this fact, since "traditionally" 
until now all the solutions obtained for both "constants" (when they vary simultaneously) 
are decreasing on time t. In such works it was not considered the possible effects of a c- 
variable into the curvature tensor in such a way that the traditional FRW equations remain 
unalterable. 

In the present paper we try to generalize the study begun in ([1]) considering a cosmolog- 
ical model with symmetries LRS Bianchi I with variable constants and which momentum- 
energy tensor is modeled by a bulk viscous fluid (full causal theory) in such a way and 
expecting that the thermodynamics help us to restrict the parameters to find physical solu- 
tions. Furthermore bulk viscosity is expected to play an important role in the early evolution 
of the Universe, when also the dynamics of the "constants" G, c and A could be different. 
We shall consider the possible effects of a c-variable into the curvature tensor in such a way 
that we shall outlined a new field equations and we shall study the Krestchmann invariants 
as well as the expansion and the shear. Once the field equations have been outlined we will 
need to make simplifying hypotheses in order to try to obtain a complete solution to the 
field equations. 

Under the assumed hypotheses we will study three different models finding that the 
physically reasonable models (decreasing temperature, increasing entropy, negative viscous 
pressure etc..) are those with G and c as growing function on time t and in particular such 
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models have an equation of the state lo — 1 that is to say, ultrastiff matter, and the viscous 
parameter is 7 = 1/2, the usual one. In such solutions the cosmological "constant" is a 
negative decreasing function. The models which have decreasing functions G and c result 
without any physical sense i.e. positive viscous pressure, decreasing entropy etc.... We end 
the paper summarizing all these results. 



II. THE MODEL AND THE HYPOTHESES 



A locally rotationally symmetric (LRS) Bianchi I ([2, 7]) spacetime with metric 

ds 2 = -c 2 (t)dt 2 + X 2 (t)dx 2 + Y 2 (t) (dy 2 + dz 2 ) , (1) 

filled with a bulk viscous cosmological fluid with the following energy-momentum tensor 
[8, 9]: 

Tt = (p + p + U)u iU k -(p + U)St (2) 

where p is the energy density, p the thermodynamic pressure, II is the bulk viscous pressure 
and Ui is the four velocity satisfying the condition U{U l = — 1. The number 4— flux and the 
entropy 4— flux take the form: 

N* = nu\ (3) 
S l = sN* - (^] u\ (4) 



2£T \ 

where n is the number density, s the specific entropy, T > the temperature, £ the bulk 
viscosity coefficient and r > the relaxation coefficient for transient bulk viscous effect (i.e. 
the relaxation time). 

The fundamental thermodynamic tensors (2-4) are subject to the dynamical laws of 
energy-momentum conservation, number conservation and the Gibb's equation: 

A k = 0, (5) 
N% = 0, (6) 

Tds = d (pi + pd ( \ j . (7) 
The equations (2-4) and (5-7) imply: 



n/ \n 



T4 = -n(3tf + ^ri + ^Tn(^) |, (8) 



(where 3H = u\) by equations (5-7) and (8), the simplest way (linear in II) to satisfy the 
i^-theorem (i.e. for the entropy production to be non-negative, S t i = ^ > ) leads to the 
causal evolution equation for bulk viscosity given by [8] 

rli + IT = -3CH -irn^H + ^-i-^y (9) 

In eq.(9), e = gives the truncated theory (the truncated theory implies a drastic 
condition on the temperature), while e = 1 gives the full theory. The non-causal theory has 
r = 0. 

The growth of the total commoving entropy £ = snV over a proper time interval (to,t) 
is given by [8]: 

£(*)-£ (to) = -t^- f^dt, (10) 

where is the Boltzmann's constant, and V = (XY 2 ) is the proper volume. 
The Einstein gravitational field equations with variable G, c and A are: 

R lk - \g* k R = ^!r^ T ik + A(t)# fc . (n) 

Applying the covariance divergence to the second member of equation (11) we get: 

div (^T> + S~h\ = 0, (12) 

% = _ - d3) 

that simplifies to: 



G 1 8ttG 



Ac 4 G c 

p + 3 (p + p) H + 3HU = -— - p- - 4p- , (14) 
Therefore, our model (with LRS BI symmetries) is described by the following equations: 



XY Y 2 8%G 
V 

Y Y 2 Yd 8ttG 



r XY + Y- 2= ^ P + Ac ^ (15) 



2 y + h ~ r Y~c = {p + n) + Ac2 ' (16) 

? (p + n) + Ac 2 , (17) 
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r + 
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-P + n) 










f 2- 



: c 2 
Y\ Ac 4 G , c 

Fj=-8^-' ) G +4 ' , ? (18) 

- + n + f ,| + 4)=-|,n((| + 4) + i-f-|,. (19 , 
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In order to close the system of equations (15-19) we have to give the equation of state 
for p and specify T, £ and r. As usual, we assume the following phenomenological (ad hoc) 
laws [8]: 

P = up, (20) 
£ = k lP \ (21) 
T = D SP S , (22) 

r = ip- l = ky-\ (23) 

where < oj < 1, and /c 7 > 0, Dg > are dimensional constants, 7 > and 5 > are 
numerical constants. Eqs. (20) are standard in cosmological models whereas the equation 
for r is a simple procedure to ensure that the speed of viscous pulses does not exceed the 
speed of light. These are without sufficient thermodynamical motivation, but, in absence 
of better alternatives, we use these equations and expect that they will at least provide an 
indication of the range of possibilities. For the temperature law, we take, T = D$p 5 , which 
is the simplest law guaranteeing positive heat capacity. 

In the context of irreversible thermodynamics p, p, T and the particle number density n 
are equilibrium magnitudes which are related by equations of state of the form p = p(T, n) 
and p = p(T,n). From the requirement that the entropy is a state function, we obtain the 
equation 



dp\ p + p T f dp 



dn ) T n n \ dT 



(24) 



For the equations of state (20-23) this relation imposes the constraint 5 = ^-j- so that 
< 5 < l/2for0<u;< 1, a range of values which is usually considered in the physical 
literature . 

The Israel-Stewart-Hiscock theory is derived under the assumption that the thermody- 
namical state of the fluid is close to equilibrium, that is, the non-equilibrium bulk viscous 
pressure should be small when compared to the local equilibrium pressure |IT| << p = up 
. If this condition is violated then one is effectively assuming that the linear theory holds 
also in the nonlinear regime far from equilibrium. For a fluid description of the matter, this 
condition should be satisfied. 

Therefore, with all these assumptions and taking into account the conservation principle, 
i.e., div(T-) = 0, the resulting field equations are as follows: 



XY Y 2 

2 1 

XY Y 2 



87TG 



p + Ac 2 




(25) 
(26) 
(27) 

(28) 

(29) 
(30) 



where W = 

Since we have defined the 4— velocity as: 

1 



u 



0,0,0 
c(t) 



U Un 



(31) 



then the expansion is defined as follows 



u. 



and the shear is 



1 



cr 



:Oi A (7 



a 




(32) 



(33) 



Curvature is described by the tensor field R l j H - It is well know that if one uses the singular 
behaviour of the tensor components or its derivates as a criterion for singularities, one gets 
into trouble since the singular behaviour of the coordinates or the tetrad basis rather than 
the curvature tensor. To avoid this problem, one should examine the scalars formed out of 
the curvature. The invariants K\ and Ki (the Kretschmann scalars) are very useful for the 
study of the singular behaviour: 



RumR 1 "* — — 7 




YYc 



(34) 
(35) 
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K 2 : = RijBP = - 



xY ( cxY X 2 cX (XY^' 

x) + \7x) - 2 x7x +3 \xy 




(Yc\ 2 YYc XY XY c YcX XYfc\ 2 

+ 3 [Ye ) - 6 YYc + 2 XY - 2 XYc ~ 2 Y7x + 2 XY [c ) + < 3?) 

+2 1L1 - 2 (*)'Z.i + 2 r(t)' + - 2* ( iV f - (38) 

XY X \Xj Yc Y\Y YYX X \Y J c K ' 




(39) 



A. Simplifying hypotheses 

In order to solve these equations we will need to make the following simplifying hypothe- 
ses: 

1. We assume a relationship between the spatial directions 

X w Y n (40) 

with n^l, and we will impose that the bulk viscous pressure follows a similar behavior 
as the energy density i.e. 

n « p, (41) 

since [II] = [p] . 

The bulk viscous pressure and the energy density of the cosmological fluid are propor- 
tional to each other and hence, the general evolution of II is qualitatively similar to the 
evolution of the thermodynamic pressure p, both obeying a similar equation of state. 
In fact, by defining an effective coefficient u e ff = u — x, the equation of state of the 
cosmological fluid can be formally written as p e ff = p + II = (u e ff — 1) P- However, 
since the bulk viscous pressure of the cosmological fluid must also satisfy the evolution 
equation (30), the resulting time evolution depends not only on uj e ff, but also, via the 
coefficients S and 7, on the equations of state of the bulk viscosity coefficient, £ = £ (p) 
and of the temperature, T = T (p) . Therefore, even that formally one can introduce 
an effective pressure (including both p and II), obeying a tu-law equation of state, due 
to the extra-constraints imposed by the requirements of the causal bulk evolution, the 



general dynamics of the present model cannot be reduced to a perfect fluid model 
evolution. 

2. As in previous works [2, 4], but here, with the difference of considering a c- variable, 
we see from the equations (26) and (27) that it is obtained 

Y Y 2 fx Y\ c XY X 

- + - + ^-- ¥ j------ = (42) 

(a) then we can make the next hypothesis: we impose that (as in the above case) 

X « F n , (43) 

therefore equation (42) yields 
obtaining a trivial solution 



Y = ( C 2 {n + 2) + C 1 (n + 2) / c(t)dt) (45) 

then we will need to make a hypothesis about the behaviour of c(t) in order to 
obtain a complete solution. 

(b) if we make the next assumption 

c(t) = c Y n (46) 
then eq. (42) has the following solution: 

Y = X (c 2 (n - 2) J X n ~ 3 dt + d(2 - n)J (47) 

then we will need to make a hypothesis about the behaviour of X(t) in order to 
obtain a complete solution. 

In the subsequent sections we shall study three different models which correspond to each 
one of the hypotheses. 
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III. SOLUTIONS 



A. Model 1. Hypotheses 1. 

We formalize this relationship in the following way, II = xp, and under physical consid- 
erations we take x e Mr and |x| <C 1 [13, 14]. 

In this section, we integrate the field equations with these assumptions and taking into 
account the conservation principle. This leads us to the following relationship, since 

X « Y n p + (uj + 1 + x) p(n + 2) ^pj = 0. (48) 

This trivially lead us to the well known relationship between the energy density p and the 
scale factor Y 

p = ^y-("+i+*)(n+2) Qr p = A w Y~ a , (49) 

where a = (u + 1 + x) (n + 2). It is observed that if n = 1, we recover the standard FRW 
model. 

Now, taking into account the eq. (30) and simplifying it, we obtain, 

xp (n + 2) 1 / (n + 2) \ 
k-yp^ 1 a 2 \ a J 

thereby obtaining p = p(t). On simplifying further, we obtain, 

p K_ 

where 



P = dk~H\ (51) 



K = I <n ^ . " . sr- I, rf=( 7 -l)J^and 6=—-. (52) 



(n+g) , (n+2)x j_ W 
a ' 2a 

From equation (49) we obtain: 



A \ 1/q 
— k b f-' 



Y=(^k b t~ b ) , i.e. Y oct(-+ 1 +-)^- 1 )("+ 2 ). (53) 



An important observational quantity is the deceleration parameter g = ^ (-^) — 1, where 
H = ^ J2i = | + , in this case = ^j^y. The sign of the deceleration parameter 
indicates whether the model inflates or not. The positive sign of q corresponds to "stan- 
dard" decelerating models whereas the negative sign indicates inflation. In our model, the 
deceleration parameter behaves as: 

q = -1- 77 -r. (54) 

H b(n + 2) y 1 



Thus, for the parameters 7, u and x, the deceleration parameter indicates an inflationary 
behaviour when 3a/b(n + 2) > and a non-inflationary behaviour for ajbin + 2) < —1. 
We proceed with the calculation of the other physical quantities as under: 



£ = fc 7 p 7 oc (dk- b t b Y = tTfeJ-^r 6 , 
T = D s p 5 = D s (dk~ b t b ) 5 with 5 = 



UJ + 1 



T = e/T 1 = A; 7 (dk- b t b ) h 1} , i.e. r = d^H. 



(55) 
(56) 

(57) 



We see from r = ef _1 £ that this result is in agreement with the theoretical result obtained 
in [8] . For viscous expansion to be non-thermalizing, we should have r < t, or otherwise the 
basic interaction rate for viscous effects should be sufficiently rapid to restore the equilibrium 
as the fluid expands. 

The comoving entropy is taking into account the relationship 



- S (t ) 



x(n + 2) k^^ 
d ^ +S -i ak B D s L 



£&(l-2±i-«) 



(58) 



n + 2 



(59) 



S(t)-S(t ) = S p (1 

We notice that the parameter x weakly perturbs the perfect fluid FRW Universe. When 
x = 0, the comoving entropy assumes a constant value and we recover the perfect fluid case. 
Finally, we will use the equations 



87rG A 2 

Ac 4 G c 

+ 4- = 0, 
mGp G c 



(60) 
(61) 



to obtain the behaviour of the "constants" G, c and A. From (60), we obtain A and using it 
in eq. (61) we obtain 

/o„ 1 1 \ o2 „2±—2—b Ta 11 I. 

(62) 



(2n + 1) /9 2 c 2 t 



47TG? 



Gfc- & 



c 1 



0. 



where /3 = —b/a. This is a Bernouilli ode and its solution is: 



c(tf 



2bk; 



t 2 (2n + 1) A 



G(t)t 3+b dt 



(63) 



r,2 

with A = f^. We can see that the above equation verifies the relationship 



k b c 2 

4- = gt b+2 therefore G = gk'cH 
G 



b 2j_-b-2 

7 



(64) 
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g being a numerical constant. It is observed that if b = —2 i.e. 7 = 1/2, we obtain G/c 2 = B. 
We use this relation in eq.(62) obtaining 



where K c = ^ 2n ^g , and therefore, 



where K c is an integration constant and K = 
Therefore the expansion and the shear yield 



■1 r . 



c = K c t\ 



(65) 



(66) 



4 _ JL 



and 



VS fx Y 



a 



3c(t) \X Y 
while the Kretschmann scalars are: 



b (n + 2) 

y/3 b(n- 1) 
3 K r at K+1 



K\ « Cit- 4{K+1) , if 2 « C 2 ^ 4(k+1) 



(67) 



(68) 



(69) 



we have omitted all the tedious calculations. 

If we impose a new assumption (which is supported by observational considerations of A 
decaying with time and assuming a small but non-zero value) as 



A 



I 



c 2 (t)t 2 



A 



21 ft r 

C**2 \C + t ' 



(70) 



we can see from eq. (60) that 

G 



3f3 2 -l 
87id 



6 6_ 2 = sjfe 6 r 6-2_ 



(71) 



Taking all these relations into consideration, we use in eq.(61) to obtain 

/ (c 1' 



whose trivial solution is 



c b + 2 

2 = 0, 

AitB \c t c t 



(72) 



(73) 



where & c is an integration constant and ji = y ' + t+g^~ 2 ^ ) • 
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1. Conclusions for this model. 

We have solved this model under the assumptions: div(T) = 0, II = xp, with x e M~ 
and valid for V7. It is observed that if 7 = 1/2 we obtain G/c 2 = const, as obtained as 
integration condition in our previous paper [1]. 

In the considered model the evolution of the Universe starts from a singular state, as 
the Krestchmann scalars show us, with the energy density, bulk viscosity coefficient and 
cosmological constant tending to infinity. At the initial moment t — the relaxation time 
is r = 0. Generally the behavior of the gravitational constant shows a strong dependence 
on the coefficient 7 entering in the equation of state of the bulk viscosity coefficient. For 
7 = 1/2, G/c 2 is a constant during the cosmological evolution. From the singular initial state 
the Universe starts to expand, with the scale factors X and Y as a 7-dependent functions. 
Depending on the value of the constants in the equations of state, the cosmological evolution 
is non-inflationary (q > for H < —1) or inflationary (q < for H > 1). Due to the 
proportionality between bulk viscous pressure and energy density, both inflationary and non- 
inflationary models are thermodynamically consistent, with the ratio of IT and p also much 
smaller than 1 in the inflationary era. Generally, the cosmological constant is a decreasing 
function of time. 

We also notice that the x— parameter, i.e. the causal bulk viscous effect, weakly perturbs 

the FRW perfect fluid solution in such a way that the comoving entropy varies with t, while in 

the perfect fluid case, i.e. x = 0, the comoving entropy is constant. The same considerations 

1 

could be taken into account for the scale factor Y(t) = Yo^ 3(1 " 7)( " +1+><)(ri+2) > which is perturbed 
weakly for the bulk viscous parameter as we can see in the special case of 7 = 1/2, and for the 
n— parameter. The solutions obtained for the scale factor, energy density, entropy, etc. are 
similar "but different" to those obtained for a perfect fluid model where the x— parameter 
vanishes. We would like to point out that our model is thermodynamically consistent for 
the usual matter equations of state and valid for all 7— parameter, i.e. V7 6 (0, 1) . In the 
same way, we can see that the viscous parameter helps us to get rid of the so-called entropy 
problem since in this model entropy varies with t. 
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B. Model 2. Hypotheses 2a. 



If we follow the hypothesis 2a then it is observed that 

Y = (c 2 (n + 2) + Ci(n + 2) J c(t)dt 

and making 

c(t) = c t a 

with o6l, therefore 



(n+2) 



y= (C 2 (n + 2) + C , 1 ^-^t a+1 , X= (C 2 (n + 2) + C , 1 



a+1 

defining 



x (c 2 (a + 1) + 

we can calculate the deceleration parameter q 

dt\H J \ C\ a + 1 J 

as we can see this expression only has sense if C 2 = 0, therefore 

2 - a 

5 



a + 1 

Now taking into account the expression 

p + (p + v + n) # = o, 

we take II as 

rr P 

n = -ap - — 

a = (a; + 1) , now simplify this expression into 

ii + -r^— + p h = -in (h + w' p - 

V 7 " 1 2 V P 

with = ( 2 ^+i ) 5 we obtain a second order ode for p(t). 

p= (§- 8 j~- ah )p- bh2 p + ch p 2 ~~' + D P 



2 
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where A = (± + *f- + a) , B = (l - § ) , C = ,D = \, simplifying it is obtained: 

p=(^-(A + l)H-f)p- BH 2 p + CHp s+1 + D— (84) 

(where 5 = 1 — 7) which is a second order differential equation with linear symmetries. 

Taking into account the standard Lie procedure [15, 18] to obtaining the symmetries of 
this ode we see that (84) only admits the following symmetry 

i = st, V = P (85) 

iff C2 = 0, (as we already know) obtaining in this way the following invariant solution: 

dt dp _ l ,1 . , 

' p = d t\ b = (86) 



( 7 -l)(f) P " U ' 7-1 
therefore 

n = -ap - =^ n = U t\ n = d k~ b f-a - (87) 

it is observed that we have obtained the same behaviour (in order of magnitude) than the 
energy density, showing in this way that our first hypothesis at least has mathematical 
meaning. We observe that this solution only has physical sense if a > 0, since U < 0. 
The entropy is 

S(t) - S (t ) = S p(i-*)+(a+i)]* o ( 88 ) 

if a > 0, then E is a growing function on time t. With these restriction we find that q behaves 
as: 

q = ^—^ < a > 2, (89) 

a + 1 

finding that this universe accelerates. 
Now from eq. 

(2n+l) — = -^p + Ac 2 , (90) 

and defining 

(2n+l)^ = /(f) (91) 

we obtain A 

««)-^)^ = A (92) 
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and substituting this expression into 

Ac 4 G 



4- = 0. (93) 



8nGp G c 
we end obtaining the behaviour of the "constant" G as: 

G= lL[L 2-\ (94) 

therefore it yields 

G = G t 2 ^- b (95) 
it is observed that we have again the relationship 



G G 

— sa t~ 2 - b , — = const. b = -2 <=^ 7 = 1/2. (96) 

c 2 & 



Once it has been obtained G then we back to eq. (92) obtaining in this way the behaviour 
of A : 

A = A r 2(a+1) (97) 

The behaviour of 

KiWt -4(a+i) ) K 2 ^t~ A{a+1 \ (98) 

while the expansion and the shear behaves as: 

e= a + l r(a+1) ^ (9Q) 



a 



K 

V3(n- 1 



(a+l)r (a+1) . (100) 



3 (n + 2) 

1. Conclusion for this model. Numerical values and graphics for the main quantities. 

In the first place we begin studying the behaviour of the radius of our model. For this 
purpose we need to fix some numerical constants as well as consider different equations of 
state. 

In the rest of this section we shall consider the following values for the numerical constants 
d = l, C 2 = 0, c = l, do = l, b= -2^7 = ^ 
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we have chosen the last values (black color) as pathological. 

With these numerical values we can see in fig. (1) the different behaviour of our scale 
factors. In all the studied cases we can see that these scale factors are growing functions on 
time t. These solutions are singular as the Krestchmann invariants show us. 





FIG. 1: In the figure (a) it is plotted the radius Y(t) while in figure (b) it is plotted the variation 
of the radius X(t). Both radius have a singular start and are growing functions on time. 

With regard to the energy density and the bulk viscous pressure (see fig. (2)) we can see 
that all the solutions have no physical meaning except the case a — |,n = | and u — 1 
(ultrastiff matter) which correspond to the red plot, since this is the only solution that 
verifies the condition II < for all time interval and for all values of the parameters. It is 
observed that for uj = (magenta line) II = 0, i.e. II vanishes. Since we have only been 
able to obtain a particular solution (invariant solution) to the differential equation (84) for 
the energy density then this quantity shows the same behaviour for the chosen numerical 
constants. 

The behaviour of the thermodynamical quantities is showed in fig. (3). Even though all 
the studied cases show a decreasing temperature except the case uj = 0, magenta line, which 

1(3 



rho(t) 




FIG. 2: Figure (a) shows the variation of the energy density pit). Figure (b) shows the variation 
of the bulk viscous pressure II (t). 

correspond to T(t) = const., as it is expected, figure (a) shows us that only the red solution 
has physical sense as already we know from the last picture of the viscous pressure. This 
case corresponds to a growing entropy. 



Sigma(l) 




FIG. 3: Figure (a) shows the variation of the entropy while figure (b) shows the variation of 
the temperature T(t). 

The variation of the "constants" G and c as well as the relationship G/c 2 is shown in fig. 
(4). This figures show us that both "constants" are growing functions on time t except in the 
case (black line) n = —1/2 where G and c are a decreasing functions on time t. We would 
like to emphasize that only for the case 7 = 1/2 we have that it is verified the relationship 
G/c 2 = const. It is a very plausible hypothesis that these effects were much stronger in 
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the early Universe, when dissipative effects also played an important role in the dynamics 
of the cosmological fluid. Hence the solution obtained (red line) could give an appropriate 
description of the early period of our Universe. 




FIG. 4: In the figure (a) we show the variation of "constant" G(t). In the figure (b) we show the 
variation of the "constant" c(t) and in the figure (c) it is plotted the relationship G/c 2 . 

With regard to the cosmological "constant", see fig. (5), it is observed that all the 
solutions are decreasing but negative except for the black line (n = —1/2). Maybe this fact 
tells us that we may consider A(t) as a true ghost energy density. 

The expansion and the shear behave as follows, see fig. (6). As we can see all the models 
studied show a decreasing expansion and the only model that has a positive shear is the 
plotted with the red color. The shear and expansion scalars calculated above indicate that 
the Universe is shearing and expanding with time. For the red color solution, the shear, 
which is a degree of anisotropy in the Universe, decreases monotonically with time. This 
indicates the fact that the initially anisotropic Universe gradually tends to an isotropic 
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FIG. 5: The variation of the cosmological "constant" A(t). 



Universe at late time (present epoch) which is in agreement with the recent observations of 
a negligible amount of anisotropy present in the CMBR. 



theta(t] 




FIG. 6: In figure (a) it is plotted the expansion 6(t), while in figure (b) it is plotted the shear a(t). 

Other solutions can be studied changing the numerical values of the constants as well as 
taking into account different equations of state. We have found at least a physical solution 
(red line) which equations of state are u — 1 (ultrastiff matter) and 7 = 1/2 (the standard 
value for the viscous parameter). 
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C. Model 3. Hypotheses 2b. 



If we follow the hypothesis 2b then it is observed that if we make the next assumption 

c(t) = c Y n (101) 

then eq. (42) has the following solution: 

Y = X (c 2 (n - 2) J X n ~ 3 dt + d(2 - nU (102) 

then we will need to make a hypothesis about the behaviour of X(t) in order to obtain a 
complete solution. Now if we assume that 

X = Y m (103) 

then 

Y = (Ci (m - n + 2) t + C 2 (m - n + 2)) <— ™+ 2 ) (104) 

defining 

g =(x +2 Fj = (m -n + 2)( Clt + g2 ) (105) 
we can calculate the deceleration parameter q 

d ( 3 \ 3n . 

9=^ 77 - 1= 3 r. 106 

H dt\H ) m + 2 v ; 

Now taking into account the expression 

p + (p + p + U)H = 0, (107) 

we take II as 

7 = — nn — 

H 



n = - a p-JL (108) 



a = (u + 1) , now simplify this expression into 



rt + r^r + ^ = 4 n f# + w- ] (109) 



with = ( 2 ^li ) > we obtain a second order ode for pit). 



^=-^--AH\p- BH 2 p + CHp 2 ^ + D— (110) 
H kj I p 
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where A = (\ + *f- + a) , B = (l - § ) , C = (j^ ,D = %, simplifying it is obtained: 

p = - ( —^— + AH + f)p- BH 2 p + CHp s+1 + D— (111) 



Cit + C 2 k 1 J r r p 

(where 5 = 1—7) which is a second order differential equation with linear symmetries. 

Taking into account the standard Lie procedure [15, 18] to obtaining the symmetries of 
this ode we see that (111) only admits the following symmetry 

£=-5(Cif + C 2 ), v = p (112) 

obtaining in this way the following invariant solution: 

d -P =► p = d (C lt + C 2 f (113) 



( T -l)(Cit + C 2 ) p 
therefore 

U = -ap--^-^U = U (C 1 t + C 2 ) b (114) 

where n = — d Q (a + b (l — 7^2)) • It is observed the we have obtained the same behaviour 
(in order of magnitude) than the energy density. 
The entropy is 

£(*) - £ (t ) = £ \(C 1 t + C 2 ) b{1 ' S)+{ ^ ) ] 1 (115) 

Now from 

(2n+l) — = — p + Ac 2 , (116) 

where 

(2n+l)^ = /(t) (117) 

we obtain A 

^~^)w) =A (118) 

and substituting this expression into 

Ac 4 G c , . 

+ 4- = 0. 119 
87rGp G c 

we end obtaining the behaviour of the "constants" . 

c 2 f ( f c\ p 
MG~p \J~ ~c)~~p 
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(120) 



G 



c 2 f f i 



Go (Cit + c 2 y 2 ~ b+ ^^ 



yf c 

with Go > 0. It is observed that we have again the relationship 
G 



'121) 



cr 



-2-b 



G 

— = const. 



7 = 1/2. 



122) 



Once it has been obtained G then we back to eq. (118) obtaining in this way the behaviour 
of A : 

-2(m+21 

(123) 



-2(m+2) 

A = A {Cit + C 2 ) m - n+2 



with A = const. 

The behaviour of the Krestchmann scalars are: 



4(u+l) 

K x « (C^ + C 2 )"^+5 



4(ra+l) 

iT 2 « (df+C 2 )-^+5 



(124) 



while the expansion and the shear behaves as: 

e = C 1 {m + 2) {C l t + C 2 ) 



_m + 2 

m — n + 2 



(125) 



a=^j-C 1 {m-l){C 1 t + C 2 ) 



m + 2 
m — n + 2 



(126) 



i. Conclusion for this model. Numerical values and graphics for the main quantities. 

As in the above case we begin fixing the values of numerical constants and the equation 
of state. In this occasion we take for the rest of this section the following values for the 
numerical constants 



and 



Ci = 1, C 2 = 1, c = 1, d = 1, b 



n 


m 
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we have chosen the last values (black color) as pathological. 

With these numerical values we can see in fig. (7) the different behaviour of our scale 
factors. In all the studied cases we can see that these scale factors are growing functions on 
time t. These solutions are non-singular as the Krestchmann invariants show us. 




2468 10 ° 2468 10 



FIG. 7: In the figure (a) it is plotted the radius Y(t) while in figure (b) it is plotted the variation 
of the radius X(t). Both radius have a nonsingular start and are growing functions on time. 

With regard to the energy density and the viscous pressure (see fig. (8)) we can see 
that the only solutions with physical meaning are those which have been plotted with red 
and blue colors i.e. with n = |,m = | and u = 1 (ultrastiff matter) for the red color 
and n — |,m = | and u — | for the blue line. Since we have only been able to obtain 
a particular solution (invariant solution) to the differential equation (111) for the energy 
density then this quantity shows the same behaviour for the chosen numerical constants. 

The behaviour of the thermodynamical quantities is showed in fig. (9). Even though 
all the studied cases shown a decreasing temperature except the case cu = 0, magenta line, 
which correspond to T(t) = const., as it is expected. Picture (a) shows us that only the 
red and blue solutions have physical sense as already we know from the last picture of the 
viscous pressure the other solutions show a decreasing entropy £(t). 

The variation of the "constants" G and c as well as the relationship G/c 2 is shown in fig. 
(10). These pictures show us that both "constants" are growing functions on time t except 
in the case (black line) n = —1/2 which correspond to a decreasing function. We would 
like to emphasize that only for the case 7 = 1/2 we have that it is verified the relationship 
G/c 2 = const. 
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FIG. 8: Figure (a) shows the variation of the energy density pit). Figure (b) shows the variation 
of the bulk viscous pressure II (t). 

M , . 2 (61 




0.5 1 1.5 2 2.5 3 1 2 3 4 5 

t t 



FIG. 9: Figure (a) shows the variation of the entropy while figure (b) shows the variation of 
the temperature Tit). 

With regard to the cosmo logical "constant", see fig. (11), it is observed that all the 
solutions are decreasing but negative except for the black line (n = —1/2). In these cases 
and because our models are non-singular A — > ±const when t — > 0. 

The expansion and the shear behave as follows, see fig. (12). As we can see all the 
models studied show a decreasing expansion and only the models that have a positive shear 
are plotted with the red and blue colors. 

Other solutions can be studied changing the numerical values of the constants as well as 
taking into account different equations of state. We have found at least two non-singular 
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FIG. 10: In the figure (a) we show the variation of "constant" G(t). In the figure (b) we show the 
variation of the "constant" c(t) and in the figure (c) it is plotted the relationship G/c 2 . 
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FIG. 11: The variation of the cosmological "constant" A(£). 

physical solution (red and blue lines) which equations of state are to — 1 (ultrastiff matter) 
and uj = 1/3 with 7 = 1/2 (the standard value for the viscous parameter). 
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theta(t) 




FIG. 12: In figure (a) it is plotted the expansion 0(t), while in figure (b) it is plotted the shear 
a(t). 

IV. CONCLUSIONS AND SUMMARY. 

In this paper we have studied three different cosmological solutions for the presented 
model. In the first of them we have needed to make two hypotheses in order to try to obtain 
a complete solution of the field equations. The first assumption, X rs Y n , is standard 
in this class of models while the second assumption, II « p,formalizes through the next 
equality, II = xp, i.e. that the bulk viscous pressure has the same order of magnitude 
of the energy density, is more difficult of digesting. But as we have seen the other two 
solutions we have founded precisely this solution. We have to point out that, as we have 
seen, this solution is not the most general solution for the outlined differential equation. 
The relationship founded between the bulk viscous pressure and the energy density is the 
invariant solution, a particular solution. Maybe the most general solution to this ode brings 
us to obtain other relationship between these quantities but as it has been showed in ([19]) 
for the classic FRW model the general solution to the outlined differential equation in the 
latter (see the appendix of ([19])) has no physical meaning while the invariant solution (the 
classical solution) at least has physical sense. For all these reason together with the reasons 
exposed above, we believe that this hypothesis is correct. 

For the first of our models we find that this is singular and thermodynamically correct 
since the entropy is a growing function on time t, if and only if x ^ 0, while the temperature 
decreases. With regard to the behaviour of the "constants" , it is founded that the constants 
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G and c for 7 = 1/2 must verify the relationship G/c 2 = const, in spite of the fact that both 
constants vary but in such a way that this quotient remains constant for all t. 

With regard to the second of our models we can see that the red color solution which 
corresponds to uj = 1 and 7 = 1/2 is the only physical solution. This solution has been 
obtained under the assumptions that X w Y n and c ~ t a , i.e. the "constant" c(t), the speed 
of light, follows an power law dependence of time t. Both hypotheses seem reasonable. The 
red solution is singular and both "constants" G and c are growing functions on time t under 
the imposed restrictions in order to find a growing entropy E and a bulk viscous pressure 
verifying the condition II < 0. As we have commented above for this case it is founded that 
II « p, i.e. both quantities have the same order of magnitude. For the studied viscous 
parameter 7 = 1/2, and only for it, we have showed that relationship G/c 2 = const, is 
verified for both constants. The cosmological constant is a negative decreasing function on 
time t. 

The third of our solutions is non-singular and have been obtained under the assumptions 
c « Y n and X « Y m . In this occasion we have obtained two physical solutions which 
correspond to the red and blue colors, with u — 1 and u = 1/3 respectively and in both 
cases the viscous parameter is 7 = 1/2, the usual one. The behavior obtained for the main 
quantities is similar to the obtained in for the second of our models except that in this 
case the solutions are non-singular i.e. the constants have a growing behavior while the 
cosmological constant is a negative decreasing function. 

Acknowledgement I wish to acknowledge to Javier Aceves his translation into English 
of this paper. 
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